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Abstract 

Computational studies show the generation of large coexisting metallic and insulating 
clusters with equal electronic density in models for manganites. The clusters are induced 
by disorder on exchange and hopping amplitudes near first-order transitions of the non- 
disordered strongly coupled system. The random-field Ising model is used to explain the 
qualitative aspects of our results. Percolative characteristics are natural in this context. 
Our results explain the recently experimentally discovered micrometer size inhomogeneities 
in manganites. The conclusions are general and apply to a variety of compounds. 



Considerable work is currently being focussed on the experimental and theoretical 
study of manganese oxides. This huge effort was triggered by the discovery of Colossal 
Magnetoresistance (CMR) in manganites (1), where the resistivity changes by several or- 
ders of magnitude upon the application of modest fields of a few Tesla at the carrier densities 
and temperatures where metallic and insulating phases are in competition. Theoretical in- 
vestigations of simple manganite models based upon numerical simulations and mean-field 
approximations have reproduced some of the complex spin-, orbital- and charge-ordered 
phases observed experimentally. In particular, the charge-ordered (CO) CE-state of half- 
doped manganites has been recently stabilized in Monte Carlo (MC) simulations of the 
two-orbital model with Jahn- Teller phonons (2). This state is in competition with the fer- 
romagnetic (FM) and A-type states also observed in experiments at the hole density x=0.5. 
However, the curious magnetotransport properties of manganites has resisted theoretical 
understanding and a proper explanation of the CMR phenomenon is still lacking. 



Recently, important new experimental information about the microscopic properties of 
manganites has been reported. Using electronic diffraction and transport techniques, Ue- 
hara et al. have investigated the effect of Pr-doping upon the metallic ferromagnetic com- 
pound La 5 / 8 Ca 3 / 8 Mn03 (3). As the Pr density (y) increases by the replacement La^Pr, 
the system changes at y~0.35 to an insulating CO-state. At low temperature in this 
regime the unexpected coexistence of giant clusters of FM and CO phases was observed 
(3). Similar results were reported by Fath et al. using scanning tunneling spectroscopy 
applied to Lao.7Cao.3Mn03 (4). At temperatures close to the Curie temperature, thin- 
film spectroscopic images revealed a complicated pattern of interpenetrating giant metallic 
and insulating phases. The clusters found in both these experiments were as large as 0.1 
^m=1000A~250a, with a~4A the Mn-Mn distance (3,4). The metal-insulator FM-CO 
transition occurs through a percolative process among the clusters, as a function of ei- 
ther temperature or magnetic field. These results rule out the picture of homogeneously 
distributed small polarons to describe doped manganites in the CMR regime. 

The discovery of huge coexisting FM-CO clusters in a manganite single-crystal is puz- 
zling. The only theoretical framework which in principle could be used to address this issue 
is the phase separation (PS) scenario where mixed-phase characteristics, involving phases 
with different electronic densities, are natural (5). The PS ideas are indeed successful in 
describing manganites at, e.g., small hole density, where nanometer size inhomogeneities 
have been widely discussed (5), and at high densities x~l based on recent magnetic and 
transport data (6). However, the micrometer clusters at intermediate densities found in 
Refs.(3,4) appear to require an alternative explanation since the energy cost of charged 
//m-size domains would be too large to keep the structure stable. Actually, explicit nu- 
merical calculations in one dimensional (ID) models have shown that the large clusters 
in PS regimes break down into smaller pieces of a few lattice spacings in size upon the 
introduction of a nearest-neighbor charge repulsion (7). In addition, this repulsion tends 
to arrange the charge in an ordered pattern (7) -charge-density- waves or stripes- contrary 
to the random location and shape of the clusters observed experimentally (3,4). A novel 
framework involving large clusters with equal-density phases is needed to rationalize the 
results of Refs.(3,4). 
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In the absence of theoretical proposals to explain the giant FM-CO clusters it is nec- 
essary to reconsider some of the properties of the models studied thus far. Of particular 
relevance is the assumption of translationally invariant interactions, which is at odds with 
the chaotic looking appearance of microclusters in experiments (3,4). For this reason, here 
we report a computational study of manganite models which combines (i) strong coupling 
interactions, necessary to reproduce the rich variety of ordered phases of these materials, 
and (ii) quenched disorder. The latter is caused by the random chemical replacement of 
ions, such as La and Pr, with different ionic sizes. This replacement affects the hopping 
amplitudes of e g -electrons due to the buckling of the Mn-O-Mn bonds near Pr (8). Re- 
cent calculations showed that the concomitant modification of the exchange coupling Jaf 
among the t2 g -spins is likely equally important in establishing the properties of manganite 
models (2). Thus, by considering randomly chosen hopping and exchange couplings fluctu- 
ating about the non-disordered values of interest, the physics of doped manganites will be 
more properly captured. Following this procedure, here we report the natural appearance 
of coexisting giant clusters of equal-density FM and AF phases in realistic models. The 
conclusions are general and similar cluster formation is expected for a variety of compounds. 

To present our main results first consider the two-orbital model, described extensively 
in previous work (5,9). It contains (i) an electronic hopping term, regulated by amplitudes 
tab' w ith a,b=l,2 labeling the d x 2_ y 2 and d3 Z 2_ r 2 orbitals, a=x,y,z being the axes directions, 
and t^=t the energy scale, (ii) a strong FM coupling between the localized t2 g - and mobile 
e g -fermions, regulated by Jh, (hi) a direct antiferromagnetic (AF) exchange among the 
localized spins with strength Jaf, and (iv) an electron-phonon coupling between the Q2 
and Q3 Jahn- Teller modes and the mobile electrons, with strength A. The phase diagram of 
the non-disordered model was studied by standard MC simulations using classical localized 
spins and phonons (2,9). Similar results were obtained with mean-field approximations 
including Coulombic repulsions (2). The generality and rationalization of the numerical 
data described below suggest that the main conclusions are actually independent of the 
detailed properties of the competing states. Whether the phases are generated by phononic, 
magnetic, or Coulombic interactions appears unimportant. 

The focus of our studies will be on first-order transitions, which in the two-orbital 
model occur in several locations in parameter space in any dimension of interest (2,5,9). 
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However, for the disorder-induced cluster formation described below it is more convenient 
to analyze ID systems first since the two-orbital 2D lattices that can be studied compu- 
tationally are not sufficiently large. Among the possible ID first-order transitions, results 
are here reported for the transition occurring between FM and AF states at fixed x=0.5 
and large A, as a function of Jaf /t. The AF phase studied has a four-spin unit cell till) 
and a concomitant peak in the spin structure factor S(q) at q=7r/2 (2,10). The state is 
insulating, as demonstrated by the absence of Drude weight and the vanishing density of 
states at the Fermi energy. Nearest-neighbors correlations among the t2 g -spins are used to 
distinguish among the FM and AF phases. In Fig.lA the energy per site (E) vs Jaf A for 
the non-disordered model is shown. The dE/d(JAF/t) discontinuity indicates the first-order 
character of the transition at JafAIc~0-21- Disorder is introduced in t" b and Jaf such 
that Jaf / 1 becomes effectively random in the interval JafAIc-^ to JafAIc+<5- Results for 
one fixed set of couplings are shown for 5=0.01 in Fig. IB (other sets lead to similar results). 
The MC averaged correlations in Fig. IB already show one of the main results of this paper, 
namely the remarkable formation of coexisting large FM and AF clusters in the ground 
state, typically of order 10a each (a is the lattice spacing). This occurs even though Jaf / 1 
at each link (not shown) rapidly changes at the a scale since different sites are uncorrelated 
in the disorder. Naively it may have been expected that at every link either the FM or 
AF phases would be stable depending on the value of Jaf A> as it occurs for a dominant 
strong disorder. However, at weak disorder this would produce a large interface energy and 
the order parameter cannot follow the rapid oscillations of Jaf A f rom site to site. As a 
consequence, structures much larger than the lattice spacing emerge, with a size regulated 
by 5 (for instance, in Fig.lC results at 5=0.05 contain FM clusters smaller than in Fig. IB). 
The effect occurs only near first-order transitions, i.e. the same weak disorder in other 
regions does not produce important effects in the spin correlations. Qualitatively similar 
results appear also in other first-order transitions of the two-orbital model, such as for the 
FM-CO(CE-state) level crossing reported in Ref.(2) using 4x4 and 4x4x2 clusters. The 
generation of large equal-density clusters by (t, JAF)-disorder near first-order transitions is 
an effect unforeseen in previous manganite investigations. 

The rapid CPU time growth with cluster size of the two-orbital model does not allow 
us to investigate numerically the phenomenon in more detail than shown in Figs.lA-C. 
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Fortunately, there are simpler models with the same behavior, including the well-known 
one-orbital model (5). It contains hopping for only one species of e g -electrons (regulated by 
t), a FM Hund coupling Jh linking the e g - and (classical) t2 g -spins, and a direct exchange 
Jaf among the t2 g -spins. Previous work showed that this model also has a first-order 
transition at x=0.5 as JafA varies, in the large Jr regime (10). It involves equal-density 
metallic FM and insulating AF states, the latter with a similar spin structure as the AF 
state of Figs.lA-C. To investigate disorder effects here the natural modification is to select 
the exchange Jaf randomly in the interval [Jaf _( Maf+^] (H)> where J AF ~0.14 is the 
critical first-order transition coupling at Jh=oo, t=l, and T=l/70 in the non-disordered 
limit (Fig. ID). In Fig. IE, results of a MC simulation corresponding to a representative set 
of random couplings (t, Jaf) centered at 0.14 are shown. As in the two-orbital case, FM 
and AF clusters, this time as large as 20a, are easily obtained. S(q) (not shown) contains 
a double-peak structure with dominant features indicating a FM-AF mixed phase. If the 
range of possible (t, Jaf) increases, the cluster size decreases (Fig. IF). Open boundary 
conditions (OBC) were used in Figs.lD-F, and periodic boundary conditions (PBC) in 
Figs.lA-C, to show that the large cluster formation occurs independently of these details. 

To further investigate the universality of the large cluster generation phenomenon, 
note that the non-disordered one-orbital model has another prominent first-order transition 
corresponding to a discontinuity in the density (n) vs chemical potential //, for a wide range 
of couplings (5). The direct interpretation of such a result is the presence of PS between 
competing FM and AF states (5). However, in the context emphasized here the focus 
shifts from the transition properties to the effect of disorder on the (n) discontinuity itself. 
Disorder is here naturally introduced as a site-dependent chemical potential of the form 
^;(/>jni, where <t>\ is randomly selected in the interval [--gr,+^]> and ni is the electronic 
number operator at site i. Results of a standard MC simulation for a L=20 sites chain of the 
disordered one-orbital model are in Fig.2A. Averages of (n) over 100 disorder configurations 
are shown. For the values of W studied here, (n) no longer has a discontinuity. For small 
W the first-order transition is replaced by a rapid crossover, where the compressibility 
proportional to d/x/d(n) remains high, suggesting the formation of large clusters. This is 
confirmed in Fig.2B where the nearest-neighbor t2 g -spin correlations are shown for a L=60 
chain, two disorder configurations, and one (typical) MC snapshot for each. These results 
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are very similar if averages over the MC configurations are made, showing that the system 
is basically frozen into an inhomogeneous ground state with large AF and FM clusters of 
size ~10-20a, similarly as in Fig.lB,E. Once again, reducing W increases the cluster sizes, 
only limited by the size of the systems that can be studied computationally. Here it was 
also observed that the mixed-phase ground state leads to a pseudogap in the T~0 density 
of states (Fig.2C), as it occurs in non-disordered models at finite temperatures in particular 
regions of parameter space (12). An analogous pseudogap was also observed working with 
the two-orbital model. Similar results as in Figs.2A-C were also found in two dimensions 
(2D) (see, e.g., Fig. 2D) and thus the large cluster formation certainly does not depend 
on pathological properties of ID systems. Note that in the particular example studied in 
Figs.2A-D, the AF-FM regions involved have different electronic densities, complementing 
the results of Figs.lA-F with equal-density clusters. In both cases the results are illustrative 
of cluster formation induced by disorder. 

The simplicity and universality of the MC simulation results with large coexisting FM- 
AF clusters suggest that there is a general principle at work in the problem. To understand 
this effect let us briefly review the phenomenology of the random-field Ising model (RFIM) 
(13) defined by the Hamiltonian H = — Jj](ij) SiSj — Xa n iSi, where Si=±l, and the rest 
of the notation is standard. The random fields {hi} have the properties [hi] av =0 and 
[h?] av =h 2 , where h characterizes the width of the distribution, and [...] av is the average 
over the fields. In manganites the Ising variables represent the competing metallic and 
insulating states on a small region of space centered at i. The random field mimics the t ab 
and Jaf fluctuations locally favoring one state over the other. Without disorder, the Ising 
model has a first-order transition at zero magnetic field and T=0 between the two fully- 
ordered states, analogous to the AF-FM first-order transitions of non-disordered manganite 
models (2). However, at h ^ the properties of the Ising transition are drastically affected 
(13). The key arguments guiding RFIM investigations (14) can be restated for manganites. 
Working very close to a first-order transition, consider that in a region dominated by 
phase-I (either AF or FM), a phase-II bubble of radius R is created. The energy cost R d_1 
is proportional to the domain wall area, with d the spatial dimension. To stabilize the 
bubble it is necessary to induce an energy compensation originated in the (t, Jaf) disorder. 
Consider the average hopping inside the bubble using Sr=^j ti, where 1 labels bonds and ti 
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is the hopping deviation at bond 1 from its non-disordered value, the latter of which is fixed 
at the critical coupling of the first-order transition of the non-disordered model. Although 
the random hopping deviations mostly cancel inside the bubble, important fluctuations 
must be considered. In particular, the standard deviation of Sr is us R =(At)R d / 2 since 
[titi']av=(At) 2 <5n/, with (At) characterizing the width of the random hopping distribution 
about the non-disordered value. A similar expression holds for the Jaf fluctuations. Then, 
centered at any lattice site it is always possible to find a region of size R, such that at least 
the average couplings favor either phase-I or -II with a substantial strength of order R d / 2 , 
although individual random deviations ti cannot exceed a (small number) At. 

To illustrate the generation of large clusters in the RFIM, standard MC simulations 
were performed. Although similar MC studies have been discussed in the RFIM framework 
(13), the results shown here provide useful qualitative information to manganite experts. 
In Fig.3A, low temperature results are shown for one representative set {hi} individually 
taken from [— W, +W] with W=3.0, in units of J=l (W=v / 3h). The dynamical formation 
of large coexisting clusters is clear, in spite of the uncorrelated character of the random 
fields in neighboring sites. Using the same set {hi} as in Fig.3A but rescaling its intensity 
with W, Fig.3B shows that as W is reduced the typical cluster sizes rapidly grow and at 
W=1.5 clusters as large as 50a in characteristic length are possible. Fig.3C contains simu- 
lation results now on a large 500x500 lattice showing that RFIM cluster sizes can be made 
as large as those found in manganite experiments (250a) by simply adjusting W. Fig. 3D 
illustrates the influence of an external field -H ex t Y2i added to the Hamiltonian. As H ext 
grows the region most affected by the field is the surface of the spin down domains, which 
are transformed into spin up. This tends to suppress the narrowest regions of the spin 
down clusters, as highlighted with arrows in Fig.3D, providing a field-induced connection 
among spin up regions that otherwise would be disconnected. Then, intuitively, as H ext 
increases a percolative transition is to be expected. Based on the RFIM-manganite analogy, 
the picture described here predicts a similar percolative transition involving metallic and 
insulating clusters as chemical compositions, temperatures, or magnetic fields are varied 
near first-order transitions, as observed experimentally (3,4). Giant cluster generation by 
weak disorder in manganite models and in the RFIM appear related phenomena. How- 
ever, at this early stage in the calculations it is difficult to predict critical exponents for 
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the metal-insulator transition. Even for simpler spin systems such as diluted anisotropic 
antiferromagnets there is still no full agreement between RFIM theory and experiments 
(13). In addition, the manganite critical dimension may be affected by the ID character of 
the zig-zag chains that form the planar CE-state (15), and critical slowing down as in the 
RFIM can produce rounding effects that make a comparison between scaling theory and 
manganite experiments difficult. 

The ideas described here are not limited to particular manganite compounds but they 
apply to other materials where a transition with first-order characteristics occurs, either by 
varying temperatures in compounds with some source of disorder, or by explicit chemical 
substitution which leads to quenched fluctuations in the hopping and exchange amplitudes. 
For instance, other manganites such as (Lai_ x Tbx)2/3Cai/ 3 Mn03 and La2-2xSri + 2 X Mn207 
also have a AF-FM competition at low temperatures. While previous investigations as- 
signed spin-glass (16) or canted-phase (17) characteristics to the intermediate region, mixed- 
phase properties involving equal-density large clusters as found in La 5 /8_ y Pr y Ca3/ 8 Mn03 
provide an alternative description. In La .5Ca .5MnO 3 , large FM and CO clusters have also 
been reported (18), and the influence of disorder on the first-order transition is a possible 
explanation for their existence (19). The concepts discussed here also apply to ruthenates, 
such as (Sri_ y Ca y )3Ru207, where a difficult-to-characterize y~0.5 region separates FM and 
AF phases (20). The metal-insulator transition of LaNii_ x Fe x 03 (21) may also proceed 
through a mixed-phase (equal-density) regime with giant cluster formation. In addition, 
EuB 6 behaves similarly to manganites (22) and it may present an analogous percolative 
behavior caused by disorder. The same could occur for the transition metal chalcogenide 
NiS2- x Se x at x~0.5 (23). Cr alloys such as Cri_ x Fe x , may also have an interesting AF-FM 
competition with percolative properties. Finally, the notorious inhomogeneities observed 
experimentally in high temperature superconductors, such as La2_ x Sr x Cu04, in regimes of 
low density of carriers may be caused in part by disordering effects on first-order transitions. 

Summarizing, based on the calculations reported here, the giant clusters in manganites 
found experimentally in Refs.(3,4) are conjectured to be caused by quenched disorder in 
the couplings (t° b and Jaf) of the system, which are induced by chemical substitution, 
and which affect transitions that otherwise would be of first-order without disorder (24). 
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A pseudogap in the density of states appears in this regime. The mixed-phase state pre- 
sented in Figs.lA-F involves clusters with equal electronic density, complementing the phase 
separation scenario which involves regions with different densities (5). In phase-separated 
regimes disorder also leads to cluster formation (Figs.2A-D). Although non-disordered mod- 
els remain crucial to determining the competing tendencies in manganites and to establish 
the order of the phase transitions, disordering effects appear necessary to reproduce the 
subtle percolative nature of the metal-insulator transition and the conspicuous presence of 
/im domains in these compounds in regimes near first-order transitions (3,4). The present 
observations are general, not based on fine tuning of models or parameters, and they should 
apply to a variety of other compounds as well. The formation of coexisting giant clusters 
when two states are in competition through first-order transitions should be a phenomenon 
frequently present in transition-metal-oxides and related compounds. 
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Figure Captions 

1. (A-C) are MC results for the two-orbital model with (n)=0.5, T=l/100, J H =oo, A=1.2, 
t=l, PBC, and L=20 (chain size). (A) is the energy per site vs Jaf / 1 for the non- 
disordered model. A level crossing (first-order transition) between FM and AF states 
occurs; (B) MC averaged nearest-neighbor t 2g -spins correlations vs position along the 
chain for one set of random t" b and Jaf couplings such that Jaf A at every site 
lies between 0.21-5 and 0.21+5, with 5=0.01. FM and AF regions are highlighted; (C) 
Same as (B) but using 5=0.05; (D-F) are results for the one-orbital model with (n)=0.5, 
T=l/70, Jh=oo, t=l, OBC, and L=64. (D) is the energy per site vs Jaf for the non- 
disordered model, showing the level crossing between FM and AF states at Jaf~0.14; 
(E) are the MC averaged nearest-neighbor t2 g -spin correlations vs position for one 
distribution of random hoppings and t2 g exchanges such that Jaf /t is now distributed 
between 0.14-5 and 0.14+5, with 5=0.01; (F) Same as (E) but with 5=0.03. 

2. Results of a MC simulation of the one-orbital model with a random chemical potential, 
PBC, J H =8.0, and J A f=0.0, in units of t=l. (A) (n) vs n for a L=20 chain at T=l/75 
using 24,000 MC sweeps per {<^} set. The results are averages over ~100 of {(f>i} 
configurations for the values of W shown; (B) Nearest-neighbors t2 g -spin correlations 
vs their location along a L=60 chain with /x=-6.7 and T=l/75. Shown are results 
for one representative MC snapshot, W=0.25 (upper panel) and W=1.0 (lower panel). 
Other snapshots differ from this one only by small fluctuations. The FM-AF clusters 
remain pinned at the same locations as the simulation evolves; (C) Density of states 
at T=l/75, L=20 and /i=-6.7 showing the presence of a pseudogap. The average 
density is (n)~0.87; (D) Results of a representative MC snapshot for an 8x8 cluster, 
T=l/50, //=-6.2 (close to the critical value), and W=1.0. Regions with FM or AF 
nearest-neighbor t2 g -spin correlations are shown. 
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3. Results of a MC simulation of the 2D RFIM at T=0.4 (J=l), with PBC. The dark 
(white) small squares represent spins up (down). At T=0.4 the thermal fluctuations 
appear negligible and the results shown are those of the lowest energy configuration. 
(A) was obtained for W=3, H ex t=0 using a 100x100 cluster and one set of random fields 
{hi}. Typical cluster sizes are ~10a; (B) Same as (A) but with W=1.5. The cluster 
sizes have grown to ~50a; (C) Results using a 500x500 cluster with W=1.2, H ext =0 
and for one configuration of random fields. The giant and percolative-like features of 
the clusters are apparent in the figure; (D) Same as (C) but now contrasting results 
between zero and nonzero H ext . The dark regions are spins up in the H ext =0 case, 
the grey regions are spins down at zero field that have flipped to up at H ext =0.16, 
while the white regions have spins down with and without the field. Special places are 
arrow marked where narrow spin down regions have flipped linking spin up domains. 
In (A)-(B) and (C)-(D) the same set {hi} was used. 
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